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Searches for high mass dilepton resonances in pp collisions at yfs = 7 
TeV with the ATLAS Experiment 

D. Olivito, on behalf of the ATLAS collaboration 

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA, USA 



The ATLAS detector has been used to search for high mass e+e - or /i + fi~ resonances, such as new heavy neutral 
gauge bosons. Over 1 fb" 1 of proton-proton collisions at ^/s = 7 TeV recorded by the ATLAS experiment at 
the Large Hadron Collider are used to search for a high mass state decaying to dilepton pairs. No excess over 
Standard Model expectations is observed, and limits are placed on benchmark models predicting spin-1 and 
spin-2 resonances. 



o 



^ 1. Introduction 

(N : 

Oh| High mass resonances are a natural place to look for new physics, and leptonic decays of these resonances would 
produce striking experimental signatures. The ATLAS detector has been used to search for high mass resonances 
decaying to leptons, where lepton in this note refers to an electron or muon. This anal ysis uses 1.08 fb -1 for 
. the e + e~ channel and 1.21 fb _1 for the channel, recorded in the first half of 2011 ( ATLAS Collaboration! 

| ' [2011aj h 

X 

i ' 

£h! 2. Theoretical Motivation 

New high mass resonances can arise in a number of models. In this note, limits are set on spin-1 and spin-2 
resonances. 

The benchmark model for spin-1 resonances is the Sequential Standard Model (SSM) Z' boson ( Langackerl 



HH), 

which shares the same production as the S tandard Model Z, as well as the same leptonic couplings 



Also considered are Eg grand unified theory models ( London and Rosnerl [198 6j | ). where the Eq group is broken 



down to SU(5) and two U(l) groups. The neutral gauge fields of the the U(l) groups, ip an d X> m i x to form 
the Z' candidate. The choice of mixing angle (0e 6 , defined in Eqn. [TJ determines the Z' coupling to fermions. 



Z'(Q E% ) = ZL cos 6 E6 + Z' sin 9 Ee ( 1 ) 



The benchmark model for spin-2 resonances is the Randall-Sundrum (RS) model ( Randall and Sundruml 



[l999j | ). which predicts excited Kaluza-Klcin modes of the graviton that appear as spin-2 resonances. These 
modes have a narrow intrinsic width when k/Mpi < 0.1, where k is the spacetime curvature in the extra 
^ . dimension, and Mpi — Mpi/\f&n is the reduced Planck scale. 



3. ATLAS Detector 



The ATLAS detector is composed of an inner detector immersed in a solenoidal magnetic field at the innermost 
radius, electromagnetic and hadronic c alorimeters, and a muqn spec trometer immersed in a toroid magnetic 
field furthest from the interaction point ( ATLAS Collaborationl |2008 |). The inner detector provides tracking of 
charged particles for pseudorapidity |»;| < 2.5. It consists of a silicon pixel detector closest to the beam pipe, a 
silicon strip detector, and a hybrid transition radiation detector and straw tube tracker at the outermost radius. 
The calorimeters span \q\ < 4.9, but for this analysis, the finely-segmented central region covering \r/\ < 2.47 is 
used. A liquid-argon electromagnetic calorimeter provides precision measurements for electrons in this region. 
Outside the calorimeters, air-core toroids provide the magnetic field for the muon spectrometer. Three sets 
of precision chambers (Monitored Drift Tubes in the barrel region, Cathode Strip Chambers in the endcaps) 
provide coverage for \rj\ < 2.7. Trigger chambers provide coverage for \rj\ < 2.4, with resistive-plate chambers 
in the barrel and thin-gap chambers in the endcaps. 
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4. Event Selection 

The goal of this analysis is to search for resonances in the high invariant mass (mi+i- ) region for electrons and 
muons. Thus the event selection requires 2 leptons of like flavor, and selection criteria are applied to the leptons 
both to minimize the reducible backgrounds and to ensure well-understood leptons at very high transverse 
momentum (pr). Both electrons and muons are triggered with single lepton triggers, with a transverse energy 
(Et) threshold of 20 GeV for electrons and pt threshold of 22 GeV for muons. For the analysis, both electron 
Et and muon pt are required to be greater than 25 GeV, ensuring maximal trigger efficiency Event candidates 
are required to have a reconstructed primary vertex with at least 3 charged particle tracks with pt > 0.4 GeV, 
to reject cosmic ray events and potential beam backgrounds. 

In addition to the above, electron candidates are required to be within \rj\ < 2.47 and are required to be 
outside the calorimeter transition region of 1.37 < |7j| < 1.52, which is not as well understood. Quality cuts 
are applied to remove electrons reconstructed in areas of the calorimeter with readout issues. As electrons 
are reconstructed from an electromagnetic calorimeter cluster matched to an inner detector track, quality 
requirements are made on both the transverse shower profile in the calorimeter and on the track. Further 
requirements are made on the matching between the two. These correspond to the Medium level of ATLAS 
electron identification ( ATLAS Collaboration! |2011b| ) and provide rejection against background arising from 



QCD jets. A hit in the innermost pixel layer is required to reduce the contribution of converted photons. 
The leading electron is also required to be isolated in the calorimeter, with ££/r(Ai? < 0.2) < 7 GeV. The 
calorimeter cells corresponding to the reconstructed electron are removed from the isolation sum, and corrections 
are applied for residual leakage of the electron energy as well as other pp interactions in the event of interest 
(i.e. pileup). The total acceptance for aZ'-> e+e~ with mass 1.5 TeV is 65%, including detector acceptance 
and efficiencies for trigger, reconstruction, and identification. The electron energy resolution is effectively flat 
at high Et, dominated by a constant term of 1.2% in the barrel and 1.8% in the endcaps. 

Muon candidates consist of a track found in the inner detector and an independent track found in the muon 
spectrometer. These tracks are combined in a single fit to determine the momentum of the muon candidate. 
To ensure high track quality, hit requirements are made in the inner detector subsystems. Specific to these 
high pt muons, the muon candidate is required to have hits in each of three sets of precision chambers in 
the muon spectrometer (inner, middle, and outer), as the momentum resolution is best understood for these 
candidates. Muon candidates are rejected if their tracks have hits in both the barrel and endcap region of 
the muon spectrometer, to minimize the impact of residual misalignments between subsystems. These cuts 
effectively restrict the geometrical acceptance. Backgrounds are suppressed by additional requirements: to 
remove contributions from cosmic rays, the muon tracks must have a transverse impact parameter consistent 
with the primary vertex, |cio| < 0.2 mm. In the longitudinal direction along the beamline, the muon track is 
required to be within |zo| < 1 mm of the primary vertex, and the position of the primary vertex relative to 
the interaction point must be within |z(.PV)| < 200 mm. QCD jet background is suppressed by requiring each 
muon to be isolated in the inner detector, such that Epx(Ai? < 0.3)/|>t(m) < 0.05, where the sum includes 
tracks with pt > 1 GeV. Furthermore, since the muon momentum measurement is taken from the curvature 
of the track, and a muon with misreconstructed charge would have grossly mismeasured momentum, the two 
muon candidates are required to have opposite charge. With these requirements, the total acceptance for a 
Z' —> ^t + /i~ with mass 1.5 TeV is 40%, again including detector acceptance as well as efficiencies for trigger, 
reconstruction, and identification. The lower acceptance relative to the electron channel is largely a result of the 
track quality requirements imposed in the muon spectrometer, which are necessary for a good understanding 
of the momentum resolution. Work is currently ongoing to understand the remaining regions of the muon 
spectrometer sufficiently well to extend this acceptance. 



5. Backgrounds 

After the selections outlined in the previous section, the primary remaining backgrounds are from Standard 
Model sources of dileptons. The main and irreducible background for both channels is Standard Model Drell- 
Yan production (Z/7*), which yields the same final state as the considered signal models. Smaller contributions 
come from electroweak diboson production (WW, WZ, and ZZ), as well as ti pairs decaying leptonically. In 
the electron channel, QCD jets can fake a prompt electron, through either a photon conversion, semi-leptonic 
decay of a b or c quark, or a hadron misreconstructed as an electron. In the muon channel, QCD jets primarily 
give rise through the semi-leptonic decays of a b and c quarks. Muons from light mesons decaying in flight 
are negligible at high pp. For both electrons and muons, jets enter the selection through either QCD dijet 
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production or production in association with a W boson (W + jets). Cosmic ray background for the muon 
channel was studied in data and found to be negligible. 

All backgrounds except QC D dijets are taken from simulation. The Z', G* signal and ZH* processes ar e 
generated with Pythia 6.421 ([Sjostrand et all [20061 ]) using MRST2007 LO* (jSherstnev and Thornd |2008| ) 
parton distribution functions (PDFs). Interference between the Z/'y* processes and th e heavy resonances i s 
small and therefore neglected. The diboson processes are generated with Herwig 6 .510 (fCorcella et al. |200l| ) 
using MRST2007 LO* PDFs. The W + .jets background is generate d with Alpgen (iMangano et al.l l2003ll"and 
the ti background with MC@NLO 3.41 ([Frixione and Webberl [20Q2jj ) . For both, Jimmy 4.31 ([Butterworth et al.l 
1996}) is used to describe multiple parton interact i ons an d Herwig to describe the remaining underlying event 



and part on showers, using CTEQ ( Pumplin et al.l [20021 ]) PDFs. Final-state photon radiation is handled with 
PHOTOS ( Golonka and Was] 120061') . The gen erated samples are p rocessed through a full simulation of the 
ATLAS detector (lATLAS Collaboration! |2010j ) based on GEANT4 (jAgostinelli et all |2003| ). 



The Z/j* cross section is calculated at next-to-next-to -leading order (NNLO) using PHOZPR (jHamberg et al.l 
199l| ) with MSTW2008 parton distribution functions ([Martin et al.1 (20091 ] V The ratio of this cross section to 



the leading-order cross section is used to determine a mass dependent QCD K-factor which is applied to the 
results of the leading-order simulations. The same QCD K-factor is applied to the Z' signal. No QCD K-factor 
is available for G* production at ^fs = 7 TeV. Higher-o rder weak corrections (beyond the photon radiation 
included in the simulation) are calculated using HORACE ( Carloni Calame et al.l [20061 [2007]), yielding a weak 
K-factor due to virtual heavy gauge boson loops. The weak K-factor is only applied to the Drell-Yan background. 
The diboson cross section is known to next-to-leading order (NLO) with an uncertainty of 5%. The W + jets 
cross section is rescaled to the inclusive NNLO calculation, resulting in 30% uncertainty when at least one 
parton with Et > 20 GeV accompanies the W boson. The ti cross sec tion is predicted at approximate-NNLO, 
with 10% uncertainty ([Moch and Uwerl [20081 ] . lLangenfeld et al.l [2009jj V 

Figure Q] shows the electron Et and muon distributions for selected leptons, with the sum of expected 
Standard Model backgrounds. The open histograms show the expected contributions from three Z' SSM mass 
hypotheses. Good agreement is observed between data and expectation. The excess of data at high px for 
muons is understood to be a combination of muon momentum resolution and poor Pythia modeling of high 
Z + jets events. The invariant mass m e + f - is modeled well by Pythia after the mass dependent QCD K-factor 
correction, however, as will be seen in Fig. [4j 
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Figure 1: Electron Et (left) and muon px (right) for selected leptons, compared with the expected Standard Model 
background. The open histograms show the expected contributions from three Zg SM mass hypotheses. 



5.1. QCD Dijets 

Due to poor modeling and low statistics in Monte Carlo, QCD dijet backgrounds are evaluated from data 
for both channels. In the electron channel, where dijet background is more significant, three independent data- 
driven methods are used. The baseline method is called "inverted identification." In it, a QCD dijet sample is 
selected in data by requiring two electron candidates which pass looser electron identification criteria but fail 
the Medium criteria used for the analysis selection. As the statistics in this sample run out before reaching 
the high mass signal region, the m e + e - distribution from this sample is fit with an empirical function to allow 
for extrapolation. To determine the normalization, a two-component fit in m e + e - is performed in the range 
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70 < m e + e - < 200 GeV, where one component is the QCD dijet function and the other is the sum of all 
backgrounds from simulation (each weighted by cross section) . 

Two other methods provide cross checks and systematic uncertainties for this estimation. The first of these, 
the "fake rate" method, derives a fake rate from a dijet sample triggered by QCD jets. To predict the background 
in the signal region, it applies this fake rate to a dielectron selection with one electron candidate passing the 
analysis cuts and the second passing a dijet-enriched selection. The last method, "isolation fits," uses a binned 
likelihood fit in calorimeter isolation to distinguish electrons from background. Since the isolation energy of jets 
increases at high Et, this method becomes increasingly effective, as long as statistics are sufficient. Templates 
are taken from single electron data: the signal template comes from electrons from W boson decays, while 
the background template comes from a jet-enriched selection where identification cuts are reversed. As seen in 
Fig. [21 the calorimeter isolation is also fairly well modeled by simulation. The leading and subleading electron 
candidates are fit separately. These results are then combined using a system of equations, to avoid double 
counting of background events. 




Leading Electron Isolation [GeV] Subleading Electron Isolation [GeV] 



Figure 2: Calorimeter isolation for leading (left) and subleading (right) electrons. 



In the muon channel, the QCD background is much smaller but still evaluated from data. A "reverse 
isolation" method is employed: a QCD sample is selected in data by requiring two muons with 0.1 < T.p^(AR < 
0.3) /pt(m) < 1-0- The normalization for this sample is obtained from the ratio of isolated to non-isolated events 
in QCD cc and bb Monte Carlo. As shown in Fig. [3l this region is dominated by QCD background and well 
described by simulation. 




Figure 3: Normalized track isolation for muons. 
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6. Results 



Figure |4] shows the resulting m^+i- distributions for electrons and muons, compared with the sum of Standard 
Model expectations. The open histograms show the expected contributions from three Z^ SM mass hypotheses. 
No clear excess is observed over the background expectations. For these plots and the limit setting procedure, 
the sum of simulated backgrounds is normalized to data in the Z peak region of 70 < mf+£- < 110 GeV, for 
reasons explained below. Table U shows the background expectations and observed data in bins of rrif+i-. In 
the first bin, data and expectation agree by construction. The QCD dijet background falls off more steeply than 
the others and is negligible in the highest m^+g- bin for both channels. 




80100 200 500 1000 2000 80100 200 500 1000 2000 

m ee [GeV] rrv, [GeV] 



Figure 4: Invariant mass (m t + t ~) for selected pairs of electrons (left) and muons (right), compared with the expected 
Standard Model background. The open histograms show the expected contributions from three Z' SSM mass hypotheses. 



Table I: Expected and observed number of events in the dielectron (top) and dimuon (bottom) channels. The first 
bin is used to normalize the total background to the data. The errors quoted include both statistical and systematic 
uncertainties, except the error on the total background in the normalization region which is given by the square root of 
the number of observed events. The systematic uncertainties are correlated across bins and are discussed in the text. 



m e+e - [GeV] 


70-110 


110-200 


200-400 


400-800 


800-3000 


DY 


258482 ± 410 


5449 ± 180 


613 ± 26 


53.8 ± 3.1 


2.8 ± 0.1 


it 


218 ± 36 


253 ± 10 


82 ± 3 


5.4 ± 0.3 


0.1 ± 0.0 


Diboson 


368 ± 19 


85 ± 5 


29 ± 2 


3.1 ± 0.5 


0.3 ± 0.1 


W+jets 


150 ± 100 


150 ± 26 


43 ± 10 


4.6 ± 1.8 


0.2 ± 0.4 


QCD 


332 ± 59 


191 ± 75 


36 ± 29 


1.8 ± 1.4 


< 0.05 


Total 


259550 ± 510 


6128 ± 200 


803 ± 40 


68.8 ± 3.9 


3.4 ± 0.4 


Data 


259550 


6117 


808 


65 


3 



"V^- [GeV] 70-110 110-200 200-400 400-800 800-3000 



DY 


236319 ± 320 


5171 


± 


150 


483 


± 22 


40.3 


± 2.5 


2.0 


± 0.3 


it 


193 ± 21 


193 


± 


20 


63 


± 6 


4.2 


± 0.4 


0.1 


± 0.0 


Diboson 


307 ± 16 


69 


± 


5 


25 


± 2 


1.7 


± 0.5 


< 


0.05 


W+jets 


1 ± 1 


1 


± 


1 


< 


0.5 


< 


0.05 


< 


0.05 


QCD 


1 ± 1 


< 


0. 


5 


< 


0.5 


< 


0.05 


< 


0.05 



Total 236821 ± 487 5434 ± 150 571 ± 23 46.1 ± 2.6 2.1 ± 0.3 



Data 236821 5406 557 51 5 



G 
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Table II: Summary of the dominant systematic uncertainties on the expected signal and background yields at m e + e - — 
1.5 TeV for the Z' (G*) analysis. NA means not applicable. 



Source dielectrons dimuons 

signal background signal background 



Normalization 


5% 


NA 


5% 


NA 


PDFs/as 


NA 


10% 


NA 


10% 


QCD K-factor 


NA 


3% 


NA 


3% 


Weak K-factor 


NA 


4.5% 


NA 


4.5% 


Trigger/Reconstruction 


negligible 


negligible 


4.5% 


4.5% 


Total 


5% 


11% 


7% 


12% 



7. Uncertainties 

Since the expected backgrounds are normalized to the observed data in the Z peak region, luminosity and 
other uncertainties which are independent of mg+£- (e.g. efficiencies) will cancel out between the Z and Z'. 
The remaining uncertainties are primarily theoretical. The 5% uncertainty on the Z cross section is taken as 
a normalization uncertainty. At higher m e + e - , the theoretical uncertainties on the Z/j* cross section grow. 
The overall uncertainty due to PDF and as variations is estimated to be 10% at 1.5 TeV using the MSTW 
2008 eigenvector PDF sets and other PDF sets corresponding to variations of as- The difference with respect 
to CTEQ is included as an additional 3% uncertainty. The uncertainty on the QCD K-factor is 3%, evaluated 
from variations of the renormalization and factorization scales b y factors of two around the no minal values. A 
systematic uncertainty of 4.5% is attributed to EW corrections ( ATLAS Collaboration! |201ld |). 



Experimental effects on efficiency, scale, and resolution were studied in data by using leptons from Z bosons 
to compare with simulation, then simulation was used to extrapolate these effects to the high m^+i- region. For 
electrons, energy scale and resolution were studied by comparing the Z line shape in data and simulation, and 
the resulting uncertainties are negligible for this analysis. Electron efficiencies for trigger, reconstruction, and 
identification were studied using the tag and probe method, resulting in an uncertainty that changes negligibly 
with m e + e -. An additional 1.5% identification uncertainty is applied to account for potential inefficiencies at 
high mass due to detector and radiative effects on the calorimeter isolation. 

For muons, the resolution is studied by comparing the Z line shape from data with simulation, and by 
comparing muon spectrometer and inner detector tracks. Simulation was adjusted to the data observations 
from this low mass region, and the derived uncertainty has an effect of less than 1.5% on the event yield. 
Trigger, reconstruction, and identification efficiencies were studied using the tag and probe technique. The 
combined uncertainty is estimated to be 4.5% at a mass of 1.5 TeV, dominated by a conservative estimate of 
energy loss from muon bremsstrahlung in the calorimeters. 

The dominant systematic uncertainties are collected in Table [Til with uncertainties less than 3% neglected. 
No theoretical uncertainties are applied to the signal models considered. 



8. Statistical Method 

A template fit method is used to quantify the agreement of the observed data with signal and background 
expectations. Signal mass hypotheses from 0.13-2.0 TeV were considered. For each signal hypothesis, a likelihood 
function is defined as the product of the Poisson probabilities for all m e + e - bins in the search region, computing 
the probability to obtain the observed data given the background plus signal. To quantify the agreement with 
a background-only hypothesis, a p-value is computed, representing the probability of seeing an excess at least 
as significant as the largest observed in data, in the absence of signal. A scan is performed over Z' cross section 
and mass, with no significant excesses found. The electron and muon channels yield p- values of 54% and 24%, 
respectively. 



9. Limits 

As no signal is o bserved, 95% confiden ce level (C.L.) limits are set on the number of Z' (G*) events using a 
Bayesian approach ([Caldwell et al.l [2009] ) with a flat prior for the signal cross section. Systematic uncertainties 
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are incorporated as nuissance parameters and marginalized, with a total effect of ~1% on the limits set. Limits 
on the number of signal events are turned in a limit on the cross section times brancing ratio (crB) by scaling by 
aB(Z —> 11) over the observed number of Z events. The expected limits are evaluated with pseudo-experiments 
containing only background events, calculating the 95% C.L. upper limit for each pseudo-experiment. The 
median of the distribution of limits is taken as the expected limit, and the distribution is used to determine 
the 68% and 95% contours corresponding to the expected limit ±l,2o\ Figure [5] shows the expected and 
observed limits for spin-1 (left) and spin-2 resonances (right), combining the electron and muon channels. The 
predicted aB for the benchmark spin-1 model, the Z' sfiMl is shown with its theoretical uncertainty (as the width 
of the curve), along with the i?6-motivated Z' models with extremal aB. For spin-2, the RS graviton aB is 
displayed for couplings 0.01 < k/M Pl < 0.1, with the theoretical uncertainty on k/M Pl =0.1. The observed 
and expected limits for the benchmark models are given in Table IIII1 with observed limits for the additional 
models in Table ITVl 



io- 1 r 
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Figure 5: Combined limits on aB for spin-1 (left) and spin-2 (right) resonances. The curves show theoretical cross-sections 
for various models. 



Table III: Observed (Expected) 95% C.L. mass lower limits in TeV on Zssm resonance and G* graviton (with k/ M p;=0.1). 



Model e + e~ ^ + \x~ 

Zg SM 1.70 (1.70) 1.61 (1.61) 1.83 (1.83) 

G* 1.51 (1.50) 1.45 (1.44) 1.63 (1.63) 



Table IV: 95% C.L. lower limits on the masses of i?6-motivated Z' bosons and RS gravitons G* for various values of the 
coupling k/M pi. Both lepton channels are combined. 





E 6 Z' Models 


RS Graviton 


Model / Coupling 


ryl ryl ryl ryl ryl ryl 

^1/) ^-n 


0.01 0.03 0.05 0.1 


Mass limit [TeV] 


1.49 1.52 1.54 1.56 1.60 1.64 


0.71 1.03 1.33 1.63 



10. Conclusions 

Searches were performed for narrow high mass dilepton resonances with the ATLAS detector at the Large 
Hadron Collider using 1.08-1.21 fb -1 of data from 2011. No significant excesses above Standard Model expec- 
tations were found, so limits on aB were set for various models predicting spin-1 and spin-2 resonances. The 
95% C.L. observed mass limit for the ZgSM i s 1-83 TeV, while for an RS graviton G* with coupling k/Mpi=0.1, 
the observed limit is 1.63 TeV. 
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